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ABSTRACT 


The  pressure  distribution  predicted  by  the  modified  Newtonian 
theory  is  used  to  develop  equations  for  the  aerodynamic  forces, 
moments,  and  stability  derivatives  for  components  of  hypersonic  lift¬ 
ing  configurations.  In  conjunction  with  the  equations,  a  set  of  charts 
is  presented  to  enable  simple  determination  of  the  aerodynamic  char¬ 
acteristics  of  swept  cylinders,  swept  wedges,  spherical  segments,  and 
cone  frustums  at  zero  sideslip  and  angles  of  attack  from  0  to  180  deg. 
This  method  allows  evaluation  of  most  delta  wing-body  combinations 
without  the  need  for  numerical  or  graphical  integration.  As  an  example 
of  the  procedure,  the  theoretical  characteristics  of  a  blunt,  75 -deg 
swept  delta  wing  are  calculated  and  compared  with  experimental  results. 
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1.0  INTRODUCTION 

In  the  design  and  testing  of  lifting  re-entry  configurations,  there  is 
often  the  need  for  a  simple,  approximate  method  of  predicting  the  pres¬ 
sures  and  forces  acting  on  the  vehicle  at  hypersonic  speeds.  The  New¬ 
tonian  theory  has  proven  very  useful  for  this  purpose.  A  number  of 
studies  have  shown  the  accuracy  of  this  simple  theory  in  predicting  the 
pressures  and  forces  on  such  configurations  as  sharp  and  blunted  cones 
(Refs.  1  through  5),  circular  cylinders  (Refs.  6  and  7),  hemispheres 
(Ref.  7),  and  delta  wings  (Ref.  8).  Although  the  Newtonian  theory  is 
easily  applied  to  the  calculation  of  pressure  distribution,  integration  of 
the  pressure  over  the  body  surface  to  obtain  total  forces  and  moments 
can  be  difficult  and  time  consuming.  Hence,  the  theory  is  not  always 
used  to  full  advantage.  Design  charts  which  simplify  the  evaluation  of 
body  loads  have  been  developed  for  complete  and  partial  bodies  of  revo¬ 
lution  (Refs.  5  and  9  through  13),  elliptic  cones  (Refs.  5,  14,  and  15), 
delta- wing  components  (Ref.  13),  and  three-dimensional  bodies  (Ref.  16). 
The  charts  of  Refs.  5,  9,  11,  13,  14,  and  15  provide  total  loads  and 
derivatives  for  selected  bodies,  while  the  methods  of  Refs.  10,  12,  and 
16  apply  to  arbitrary  bodies  but  require  numerical  or  graphical  integration. 

The  purpose  of  the  present  report  is  to  extend  the  scope  of  the  pre¬ 
vious  design  charts  by  providing  additional  aerodynamic  characteristics 
and  an  increased  angle -of- attack  range.  To  avoid  a  requirement  of 
numerical  integration,  only  selected  configurations  consisting  of  typical 
delta-wing  and  body  components  are  considered.  Equations  are  derived 
for  the  pressure  distribution  on  each  component,  and  this  distribution  is 
then  integrated  over  the  surface  area  in  closed  form  to  obtain  total  forces 
and  moments.  Equations  and  charts  are  given  for  the  longitudinal  sta¬ 
bility  and  performance  coefficients  (Cn,  Ca,  Cm):  and  the  directional  and 
lateral  stability  derivatives  Cn^>  f°r  an  angle -of- attack  range 

of  0  to  180  deg  at  zero  sideslip. 

An  example  of  the  use  of  the  charts  is  given  in  the  appendix,  where 
the  aerodynamic  characteristics  of  a  blunt,  75 -deg  swept  delta  wing  are 
computed  and  compared  with  experimental  results. 

2.0  DEVELOPMENT  OF  EQUATIONS 

The  Newtonian  theory  has  been  discussed  in  a  number  of  references, 
and  only  a  brief  summary  will  be  given  here.  A  thorough  analysis  of  this 
theoretical  method  is  given  by  Hayes  and  Probstein  in  Ref.  17. 


Manuscript  received  January  1964. 
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Newton  calculated  the  force  on  a  body  by  assuming  that  the  impact 
of  fluid  particles  was  completely  inelastic  for  the  normal  component  of 
momentum  and  was  frictionless.  Thus,  the  normal  component  of  mo¬ 
mentum  is  converted  to  a  pressure  force  on  the  body  while  the  tangen¬ 
tial  component  remains  unchanged.  The  analysis  based  on  these 
assumptions  gives  the  surface  pressure  coefficient,  Cp ,  as 

Cp  =  2  cos2  T]  (1) 

where  v  is  the  angle  between  the  free- stream  velocity  vector  and  the 
inward  directed  unit  vector  normal  to  the  surface. 

At  high  Mach  numbers  the  disturbed  region  in  front  of  a  body  be¬ 
comes  very  limited  in  extent.  The  bow  shock  wave  has  approximately 
the  same  inclination  as  the  body  and  is  separated  from  the  body  surface 
by  a  very  thin,  practically  inviscid,  shock  layer.  With  this  flow  geom¬ 
etry,  the  normal  momentum  of  impinging  molecules  is  lost  inelastically 
and  the  tangential  component  of  momentum  is  conserved.  Hence, 

Newton's  analysis  is  realistic  for  this  type  of  flow,  and  the  validity  of 
the  analysis  increases  as  the  shock-layer  thickness  decreases.  For 
the  shock  wave  to  approach  the  inclination  of  the  body,  the  gas  dynamic 
equations  show  that  the  ratio  of  the  density  ahead  of  the  shock  to  that 
behind  the  shock  must  approach  zero.  The  equations  further  show  that 
for  the  density  ratio  to  approach  zero,  the  Mach  number  must  approach 
infinity  and  the  ratio  of  specific  heats  must  approach  unity.  If  these 
Newtonian  conditions  (Mm  oo  9  y  -»  1)  are  satisfied,  the  Newtonian  pres¬ 
sure  coefficient,  Eq.  (1),  is  identical  to  that  given  by  the  oblique  shock 
relations  for  the  pressure  immediately  behind  the  shock  wave. 

In  Newton's  analysis,  the  impinging  molecules  leave  the  body  surface 
along  an  unaccelerated  path.  However,  in  the  case  of  a  curved  body  with 
a  thin  shock  layer,  the  particles  are  constrained  in  the  shock  layer  and 
must  follow  an  accelerated  path.  Therefore,  for  a  correct  analysis 
Eq.  (1)  must  be  modified  to  allow  for  the  pressure  gradient  resulting  from 
the  centrifugal  forces  acting  on  the  particles.  This  correction  was  first 
obtained  by  Busemann  (Ref.  18),  and  the  rational  theory  including  the 
correction  has  been  called  the  Newton-Busemann  theory  in  Ref.  17.  How¬ 
ever,  despite  the  theoretical  correctness  of  the  Newton-Busemann  rela¬ 
tion,  the  simple  Newtonian  theory  has  been  found  to  agree  much  better 
with  experimental  data  (e.  g.  ,  Ref.  7),  and  the  equations  given  in  the  pres¬ 
ent  paper  have  not  been  corrected  for  centrifugal  effects. 

Equation  (1)  has  been  modified  by  a  number  of  investigators  to  pro¬ 
vide  a  better  correlation  with  experimental  data  for  several  classes  of 
bodies.  The  modified  forms  of  the  equation  have  the  general  relation, 

Cp  =  K  COS2  7]  (2) 
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where  K  is  a  multiplicative  factor  which  is  used  to  match  certain  limit¬ 
ing  conditions.  For  a  flat  plate  with  attached  shock  (i.  e.  ,  at  low  angles 
of  attack).  Love  (Ref.  19)  suggested  that  K  =  y  +  1  provides  better 
agreement  with  the  exact  oblique  shock  solution.  For  a  slender  pointed 
body  with  attached  shock,  best  agreement  with  exact  theory  is  obtained 
by  using  either  the  simple  Newtonian  value  of  K  =  2  or  the  value 

c 

suggested  in  Ref.  19  of  K  =  —  2p-°°  — — ,  where  Snose  is  the  surface  angle 

sin  Ono s e 

at  the  nose  and  CPnoge  is  the  exact  value  of  pressure  coefficient  for  this 
angle.  Lees  (Ref.  20)  suggested  that  for  a  blunt  body  with  detached 
shock  wave  the  Newtonian  theory  could  be  modified  to  match  conditions 

at  the  stagnation  point  by  letting  K  =  Cp  =  — 1 - - - —  ,  which  is  closely 

m  a  x  *Ioo 

approximated  by  K  =  y  +  3/y  +  1  for  large  Mach  numbers.  In  the  present 

derivations,  the  modified  form  of  the  Newtonian  approximation  as  given 
in  Eq.  (2)  will  be  used  with  an  arbitrary  value  of  K. 

The  angle,  v,  between  the  velocity  vector  and  the  surface  unit 
inner  normal  vector  n  is  determined  by  the  scalar  product  of  the  two 
vectors.  The  velocity  vector  (Fig.  1)  is 

(  i  cos  a  cos  /3  -f  j  sin  /3  +  k  sin  a  cos  /3)  (3) 

where  7,  "J,  and  k  are  unit  vectors  directed  along  the  X-,  Y-,  and 
Z-axes.  The  body  surface  may  be  described  by  the  equation 
F(x,  y ,  z)  =  0.  Then  the  inward  directed  unit  vector  normal  to  the  body 
surface  is  ’  * 


n  =  i  cos  (n,  x)  +  j  cos  (n,  y)  +  k  cos  (n,  z)  (4a) 

where  (n,x),  (n,  y),  and  (n,  z)  are  the  angles  between  n  and  the  positive 
X-,  Y-,  and  Z-axes,  respectively,  and  their  cosines  are  given  by 


Thus, 


cos  (n  ,  x )  =  — 


cos  (n,  y )  =  - 


cos  (n  ,  z  )  =  — 


_ dF  /  dx _ 

\Z{dF/dx?  +  (  <9F  /  dy  Y  +  (  <9F  /  <9Z  )2 

dF  /  dy 

V  (  <?F  /  <9x  )2  +  (  dF  /  dy)2  +  (  dF  /.  dz)2 

_ (9F  /  dz _ 

V(<9F/<9x)2  +  (c?F/<9y)2  +  (dF/dz)2 


(4b) 


COS  T) 


(5) 


—  cos  a  cos  jS  cos  (n,  x)  +  sin  /3  cos  (n,  y)  +  sin  a  cos  /3  cos  (n,  z)J 
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The  Newtonian  theory  predicts  pressures  only  on  surfaces  which  face 
the  flow.  For  surfaces  which  are  shielded  from  the  flow,  it  is 
assumed  that  the  surface  pressure  is  equal  to  the  free-stream  static 
pressure  and  Cp  =  0.  Therefore,  Eqs.  (1)  and  (2)  are  applicable  only 

for  COS  r/  >  0 . 


Force  and  moment  coefficient  nomenclature  utilized  in  the  deriva¬ 
tion  is  shown  in  Fig.  1.  The  coefficients  are  non-dimensionalized  by 
an  arbitrary  reference  area,  S,  and,  in  the  case  of  moment  coefficients, 
by  an  arbitrary  reference  length,  i.  The  moment  reference  point  of 
each  component  is  given  in  the  corresponding  figure.  The  coefficients 
are  obtained  by  integrating  the  Newtonian  pressure  distribution  over  the 


body  surface  area.  A,  as  indicated  in  the  following  general  equations 
where  the  moment  reference  point  is  at  the  origin  of  the  axes: 

Cn  = 

fn 

<IooS 

-  -f-  If 

b  A 

cos2  tj  cos  (n  ,  z  )  dA 

(6) 

cA  = 

fa 
q  S 

noo 

~  4  " * 

cos2  Tj  cos  (n  ,  x  )  dA 

(7) 

r  - 

My 

-  IT  K 

x  cos2  tj  cos  (n,z)dA  —  //  z  cos2  77  cos  (n, 

A 

x)dA]  (8) 

'-,in  _ 

Cy  - 

F  Y 

-F  ffA 

cos2  Tj  cos  ( n  ,  y  )  dA 

(9) 

%oS 

r  - 

Mz 

Tf  [ffA 

x  cos  tj  cos  ( n  ,  y  )  dA  —  f  f  y  cos2  77  cos  ( n  , 

A 

x )  d  A  J  ( 1 0 ) 

'“■n  — 

q<x,Si 

Ci  - 

MX 

£  K 

y  cos2  77  cos  ( n  ,  z )  dA  —  f  f  z  cos2  tj  cos  (n  ,  y 

A 

)dA]  (11) 

<1*0  S£ 

In  the  integration  over  the  surface  area,  it  is  assumed  that  Cp  =  0 
on  all  surfaces  shielded  from  the  flow  and  on  all  flat  surfaces  (except 
in  the  case  of  the  swept  wedge)  because  these  surfaces  are  usually  con¬ 
cealed  by  other  body  components.  Equations  and  charts  are  given  for 
components  having  vertical  symmetry  and  for  the  corresponding  flat- 
topped  components.  For  flat-bottomed  components,  the  loads  may  be 
determined  by  taking  the  difference  between  the  loads  acting  on  the  com¬ 
plete  component  and  those  acting  on  the  flat- topped  component  and  adding 
the  pressure  load  of  the  flat  lower  surface  to  the  normal  force  and 
pitching  moment. 
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2.1  DELTA-WING  COMPONENTS 

The  basic  components  of  a  delta  wing  are  the  nose,  leading  edge, 
and  wing.  The  loads  on  these  components  are  computed  in  the  fol¬ 
lowing  sections.  The  method  of  combining  the  components  to  give  a 
complete  wing  is  described  in  Section  2.  3. 


2.1.1  Spherical-Wedge  Nose 


The  nose  of  a  delta  wing  is  usually  a  spherical  wedge.  If  the  wing 
centerline  angle  e  (see  Fig.  F)  is  not  zero,  the  nose  is  not  exactly 
spherical,  but  the  error  in  force  coefficients  will  be  negligible  when  e 
is  small.  The  nomenclature  used  in  the  derivation  is  shown  in  Fig.  A. 


The  direction  cosines  of  the  inward  directed  unit  normal  vector,  n,  as 
obtained  from  Eq.  (4b)  or  by  analytic  geometry  are 


cos 

(n  ,  x  )  = 

—  COS  0 

cos  6 

j 

) 

cos 

(n  ,  y  )  = 

COS  (f> 

sin  0  / 

'  (12) 

cos 

(  n  ,  z  )  = 

—  sin  <f> 

j 

1 
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The  angle  q  between  the  free- stream  velocity  and  the  normal  vector  is 
obtained  from  Eq.  (5)  which  gives 


cos  q  =  cos  /3  (cos  a  cos  <^>  cos  6  +  sin  a  sin  (f>)  —  sin  /3  cos  cf>  sin  6 


(13) 


The  pressure  distribution  over  the  nose  is  given  by  Eq.  (2)  as 


Cp  =  K 


2 

COS  Tj 


The  value  of  at  which  the  surface  becomes  shielded  from  the  flow  is 
designated  as  <£0  and  is  defined  by  Cp  =  0  or  cos  -q  =  0,  thus 


tan  (f)Q  = 


sin  sin  6  —  cos  78  cos  a  cos  6 
cos  /S  sin  a 


(14) 


and  at  (3  =  0 


(14a) 


The  elemental  surface  area  is  given  by 

dA  =  R2  cos  cf>  d  (f>  d  6  (15) 

As  was  mentioned  previously,  it  is  assumed  that  Cp  =  0  on  the  flat 
base  surfaces  even  at  angles  of  attack  where  these  surfaces  are  not 
shielded  from  the  flow.  All  coefficients  and  derivatives  are  evalu¬ 
ated  at  fi  =  0. 


2. 1.1.1  Normal-Force  Coefficient 


The  normal-force  coefficient  is  given  by  Eq.  (6): 


Cr  =  ~  ~7r“  ff  cos  V  cos  (n,  z)  dA 
S  A 


Since  the  body  has  lateral  symmetry,  this  equation  may  be  integrated 
over  the  left  side  and  the  results  multiplied  by  2.  Then,  for  0  <  a  <  n. 


Cn  = 


2  K  R 


A  7T/2 

f  f  cos  q  sin  cos  cf>  d  d  6 

o  <f>  o 


(16) 


Substituting  the  value  for  cos  q  at  yS  =  0  from  Eq.  (13)  and  performing 
the  indicated  integrations  gives  the  normal-force  coefficient  as  a 
function  of  a  and  A.  Note  that  the  equation  must  be  evaluated  by  inte¬ 
grating  first  between  <f>  =  <f>Q  and  =  n/2  and  then  integrating  the  result¬ 
ing  function  between  6=0  and  6  =  A  since  <f>Q  is  a  function  of  6. 

Then, 

CN  — ^  =  si-n-a-  F cos  a  sin  A  A 


KR 


(f 


+  tan 


tan 


— j  +  tan  1  (sin  a  tanA)J  (17) 
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Equation  (17)  was  evaluated  for  a  =  0  to  180  deg  and  A  =  60  to  90  deg 
and  the  results  are  presented  in  Fig.  2a. 

2. 1.1. 2  Axial-Force  Coefficient 


The  axial-force  coefficient  is  given  by  Eq.  (7): 

C  a  =  —  — ~  / f  cos2  7]  cos  (  n  ,  x  )  dA 

b  A 

Then,  for  0  <  a  <  n, 

2  A  77/2 

0*A  —  - —  f  f,  cos  77  cos  co  cos  6  d  d>  d  0  (18) 

S  0 


Integration  of  Eq.  (18)  gives 


CA 

s 

2 

J_  1 

(sin2  a  sin  A  +  3  cos2  a  sin  A  —  cos2  a  sin3  A)  I 

(-2-  +  tan"1 

cos  A  \ 

K  R 

4  1 

- 

'  2 

tan  a  / 

+  2 


cos  a  tan 


(sin  a  tan  A)  +  sin  a  cos  a  sin  A  cos  A 


1 


(19) 


The  numerical  evaluation  of  Eq.  (19)  is  presented  in  Fig.  2b. 

2.1.1 .3  Pitching-Moment  Coefficient 

Since  the  force  on  any  element  of  surface  is  directed  toward  the 
center  of  curvature,  the  resultant  force  acts  through  the  center,  and 
the  moment  about  the  reference  point  is  zero. 

2.1 .1.4  Side-Force  Coefficient  Derivative 


The  side-force  coefficient  is  given  by  Eq.  (9): 


Cy  =  —  //  cos2  T)  cos  (n  ,  y )  dA 

S  A 


For  either  side  of  the  nose,  at  0  ^  a  <  n, 
Cy  = 


KR 


+  A  77/2 

f  f  cos  77  cos  d>  sin  e  d  0  d  e 

0  V  o</8> 


(20) 


where  the  upper  sign  corresponds  to  the  left  side  of  the  nose  and  the 
lower  sign  corresponds  to  the  right  side.  The  derivative  of  the  side- 
force  coefficient  with  respect  to  /3  is 


kr2 

s 


±A 

/  sin  6 


d 

d/3 


77/2 

/  cos  rj  cos  cf>  d  cf)  d  6 


(21) 
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Since  the  lower  limit,  <f>0,  is  a  function  of  p,  the  Liebnitz  rule  is  used 
to  obtain  the  derivative 


d 

dp 


n/2 

f  COS2 

V.(/3> 


V 


<f>  d  </> 


n  /  ft / 2 

=  -  cos2  77  (<£o)cos2-0o  +  f 

dP  <t>0  (j8) 


dp 


2 


COS 


<f>  d  <f> 


By  definition,  cos2  77  (<£o)  =  0,  and  d(f>o./dp  is  finite,  so  the  first  term  is 
zero.  Thus 


K  R2 
S 


+A  TT/z 

f  L  . 

°  0O (P) 


d  cos  rj 

~w~ 


2 

COS 


( p  sin  0  d  (f>  d  6 


(22) 


Since  <f>  and  6  are  independent  of  P,  the  derivative  at  P  =  0  is 

+A  JT/2 

=  ±  -K  ^  I 

P  =  o 

From  Eq.  (13) 


/  d  cos  77  j 

J0  J(/>o(P  =  0)  \  dp  J p  =  o 


(Cyp)p=„  =  *  ^r~  t  ■"  1  "  7/  - 1  -  cos'  4,  sin  e  a  <h  e 


(23) 


d  cos  77  J 

~dp  ~)p  =  0 


2  cos  <;!>  sin  0  (cos  a  cos  (f>  cos  0  +  sin  a  sin  <f>) 


(24) 


Substituting  Eq.  (24)  in  Eq.  (23)  shows  that  Cy^g  is  the  same  for  both 
sides  of  the  nose.  Then,  multiplying  Eq.  (23)  by  2  and  performing  the 
integration  gives 


s 

k  n2 


cos  a  sin3  A 

(—  +  tan 

cos  A\ 

_ 

\  2 

tan  a  / 

+  tan  1  (sin  a  tan  A)  —  sin  a  sin  A  cos  A  j 
The  numerical  evaluation  of  Eq.  (25)  is  presented  in  Fig.  2c. 


(25) 


2. 1.1. 5  Yawing-Moment  and  Rolling-Moment  Coefficient  Derivatives 

As  was  the  case  with  pitching  moment,  the  yawing-moment  and 
rolling- moment  coefficients  and  their  derivatives  are  zero  about  the 
reference  point. 


2.1.2  Flat-Topped  Spherical- Wedge  Nose 

The  geometry  of  the  flat- topped  spherical  wedge  is  shown  in  Fig.  B. 

The  direction  cosines,  pressure  coefficient,  and  elemental  area 
are  the  same  as  for  the  complete  spherical  wedge.  For  0  <  a  <  n/2  no 
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'jr- 


t 


|  i 


Fig.  B  Flat-Topped  Spherical-Wedge  Nose 


part  of  the  curved  surface  is  shielded  from  the  flow,  so  the  limits  of 
integration  on  ^  become  0  to  n/2.  Then,  Eqs.  (16),  (18),  and  (23)  give 
the  following  results: 


CN  -"I  =  —  ^ tt  sin  a  cos  a  sin  A  +  cos2  a  sin  A  cos  A  +  A  (1  +  sin2  a  )]  (26) 

cos’  a  (sin  A  -  si°’  A ) 


CA 


KR 

S 

KR2 


j_  r  3  n 
4  L  2 


77  •  2  •* 

+  —  sin  a  sin  A 


(27) 


+  2  cos  a  sin  a  ( A  +  sin  A 


S 

K  R2 


77  cos  a  sin 


A  +  2 


(A  - 


A 


A) 


(28) 


where  A  is  in  radians.  The  moment  coefficients  and  their  derivatives 
are  zero  for  the  indicated  reference  point.  For  a  >  n/2  the  equations 
for  the  complete  spherical  wedge  also  apply  to  the  flat-topped  spherical 
wedge. 
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The  characteristics  of  the  flat-topped  spherical  wedge  are  pre¬ 
sented  in  Fig.  3. 

2.1.3  Swept-Cylinder  Leading  Edge 


The  delta- wing  leading  edge  which  is  analyzed  in  this  section  con¬ 
sists  of  two  symmetrically  swept  circular  cylinders.  The  nomencla¬ 
ture  used  in  the  derivation  is  shown  in  Fig.  C. 


Fig.  C  Swept-Cylinder  Leading  Edge 

The  two  sides  are  treated  as  a  unit  in  the  analysis,  and  the  sweepback 
angle  A  is  taken  as  positive  in  the  equations  for  the  leading- edge  coef¬ 
ficients.  Where  the  two  sides  must  be  considered  separately,  as  in 
Eqs.  (29)  and  (30),  the  sign  convention  is  A  >  0  for  the  right  side  and 
A  <  0  for  the  left  side. 
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The  direction  cosines  of  the  inward  directed  unit  normal  vector  are 

cos  (n,  x)  =  —  cos  (f>  cos  A 
cos(n,y)  =  —  c os  cji  sinA 
cos  ( n ,  z  )  =  —  sin  0 

Then,  from  Eq.  (5) 

cos  ij  =  cos  £  (cos  a  cos  <£>  cos  A  +  sin  a  sin  cf>)  +  sin  /3  cos  cf)  sin  A  (30) 

At  jS  =  0  the  surface  becomes  shielded  from  the  flow  along  a  line  de¬ 
fined  by  <f>  =  4>0 ,  where 

(Ta“)  <31> 

The  elemental  surface  area  is 

dA  =  2LR  d<£  (32) 

The  integration  of  the  pressure  distribution  over  the  surface  has  as 
its  limit  the  geometric  angle  <£',  which  is  always  positive.  For  a  lead¬ 
ing  edge  which  is  tangent  to  the  wing  surface,  </>'  is  determined  by  the 
wing  sweepback  and  dihedral  angles  A  and  T,  as  shown  in  Fig.  D. 


A 


V 


Fig.  D  Leading  Edge  and  Wing  Geometry 
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Thus, 

and 

Therefore 


4> '  - 

-  y 

2  r 

(33a) 

y  = 

(if) 

(33b) 

r  = 

t“"1  (Zr) 

(33c) 

For  a  leading  edge  which  is  not  tangent  to  the  wing  surface,  Eq.  (33)  is 
not  valid  and  <j>'  must  be  determined  from  the  leading- edge  geometry. 

It  is  assumed  that  Cp  =  0  on  all  flat  surfaces,  and  the  coefficients  and 
their  derivatives  are  evaluated  at  /3  =  0. 


2. 1.3.1  Normal-Force  Coefficient 


The  normal-force  coefficient  is  given  by 

4>' 

Cn  =  -  2  K  L  jt  j  cos2  77  sin  cfi  d  cf>  (34) 

‘<f>0  (  »r-f) 

If  <j>'  <  |  </>0  I,  the  lower  limit  of  integration  is  <f>  =  -  <f>' .  The  limit 
changes  to  <f>  =  <£0  at  4>r0  =  -  or  at  a  -  a0  ,  where 


Evaluating  Eq.  (34)  for  0  <  a  <  «0  with  the  lower  limit 

CN  — — —■  -  =  sin  a  cos  a  cos  A  sin3  d> ' 

K  L  R  3 


and  for  a0 


Cn 


s 

KL  R 


<  a  <  ( 77 


-  aQ)  with  the  lower  limit  4>0  > 
a  —  cos  a  cos  A)  (cos  (f> '  —  cos3  cfiQ) 


(35) 


(36) 


(37) 


—  3  sin  a  (cos  <£> '  —  cos  <£0 )  +  2  sin  a  cos  a  cos  A  (sin3  cf>'  - 


<f>o ) 


where  4>0  and  <f>'  are  given  by  Eqs.  (31)  and  (33c),  respectively.  If  the 
wing  dihedral,  F ,  is  small,  the  leading  edge  will  be  closely  approxi¬ 
mated  by  a  complete  hemicylinder  (<£'  =  n/ 2)  .  Then,  from  Eqs.  (31) 
and  (37)  for  o  <  a  <  n 


Cn 


_ 

K  L  R 


A 

3 


sin  a 


|cos  a 


A  +  V  1 


A 


cos  a 


(38) 


Equation  (38)  was  evaluated  for  a  =  0  to  180  deg  and  A  =  60  to  90  deg, 
and  the  results  are  presented  in  Fig.  4a. 
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2. 1.3.2  Axial-Force  Coefficient 


The  axial-force  coefficient  is  given  by 


C  a  = 


f  COS  7]  COS  0  cos  A  d  0 

s  c f>0  (  or  —  0  ') 


Then,  for  o  <  a  <  a0 


Cy  v  T~”-n  =  CQS  ^  (sin2  a  -  cos2  a  cos2  A). sin3  0 '  +  3  cos2  a  cos2  A  sin  0 '  (40) 

K  L  R  3  L  -l 


and  for  a0  <  a  <  U  -  aQ) 


2  cos  A  /  •  2 


(sin2  a  —  cos2  a  cos2  A)  (sin3  0  '  —  sin3  0O) 


+  3  cos2  a  cos2  A  (sin  <f>'  —  sin  0O) 
—  2  sin  a  cos  a  cos  A  (cos3  tf>'  —  c( 


For  the  hemicylinder  (0'  =  n/2)  , 


S  _  4  cos  A 
KLR  3 


+  cos2  a  cos2  A  +  cos  a  cos  A  V  1  —  sin2  A, 


cos  a  ;  (42) 


The  numerical  evaluation  of  Eq.  (42)  is  presented  in  Fig.  4b. 

2. 1.3. 3  Pitchinq-Moment  Coefficient 


The  resultant  force  acts  through  the  center  of  curvature  at  a  point 
midway  between  the  ends.  Thus, 


C  m  =  C]\j 


L  sin  A 


2. 1.3. 4  Side-Force  Coefficient  Derivative 

The  side-force  coefficient  for  either  side  (with  proper  sign  con¬ 
vention  on  A )  is 

0' 

Cy  =  —  ^  ^  R  j  ,  cos2  rj  cos  0  sin  A  <1  0  (44 

S  -  0O-  {  or  —  0  ') 

Using  the  same  procedure  as  was  used  with  the  spherical-wedge  nose, 


Cyfl  =  —  -FTJL  f  I  ^  — 2.  )  COS  0  sin  A  d  0 

P  s  1  <3/3  =  0  v 
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where 

/ d  COS  7) 


3P  7/8- 


=  2  cos  (f)  sin  A  (cos  a  cos  cos  A  +  sin  a  sin  0) 


(46) 


Substituting  Eq.  (46)  in  Eq.  (45)  shows  that  Cy^  is  the  same  for  both 
sides.  Then,  multiplying  Eq.  (45)  by  2  and  performing  the  integration 
gives  for  0  <  a  <  aQ 


c  .2  k  (  .  ,  s  sin 

CYn  =  -  8  cos  a  cos  A  sin  A  1  sm  (p  -  - - - 1 

o  K  L  R  '  o' 


■jS  KLR 
and  for  aQ  <  a  <  (n  -  aQ) 


(47) 


Cy 


(8  KLR  3 


4  •  2  * 

sin  A 


£  cos  a  cos  A  (3  sin  cf>'  ~  3  sin  0O  —  sin3  <f>'  +  sin3  0O) 

—  sin  a  (cos3  <f>'  —  cos3 


(48) 


For  the  hemicylinder  (0'  =  n/2) , 


Cy 


/3  KLR 


A 


o _ _ _ a  .  sin2  a  +  2  cos2  a  cos  A 

Z  cos  a  cos  A  +  .  . == . =  . — 

V  1  —  sin  A  cos  a 


^he  numerical  evaluation  of  Eq.  (49)  is  presented  in  Fig.  4c. 
.1.3.5  Yawing-Moment  Coefficient  Derivative 

The  yawing-moment  coefficient  is  given  by  Eq.  (10): 


(49) 


Cn  = 


j^ff  x  cos2  V  cos  (n,y)dA—  f  f  y  cos2  rj  cos  (n ,  x)  d  aJ 

rom  Fig.  C,  x  -  ±  sin  A  and  y  =  ±  cos  A  ,  where  the  upper  sign 

>plies  to  the  right  side  and  the  lower  sign  applies  to  the  left  side, 
nee  x  and  y  are  not  functions  of  the  yawing- moment  coefficient  for 
ther  side  is 


Cn  =  ± 


2  J 


(Cy  sin  A  +  Ca  cos  A) 


mparing  Eq.  (39)  and  Eq.  (44)  gives 


Ca  = 

A  tan  A 


Cn  =  ±  Cy 


/  2  sin2  A  —  1  \ 
\  sin  A  / 


total  yawing-moment  coefficient  for  both  sides  i 
Cn  =  (CY,ight  +  CYlcIt)  (-'  S2£  ~ 
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and  the  derivative  is 


C"/S  (CYiSright  +  2  1  (' 


2  sin  A  —  1 
sin  A 


But,  CYn  =  CYn  =  — — 
P  right  P  left  2 


'/8 


Therefore 


S  -  CyS  ~TJ 


(  2  sin2  A  —  1  \ 
\  sin  A  / 


(50) 


2. 1.3.6  Rolling-Moment  Coefficient  Derivative 

The  rolling-moment  coefficient  is  given  by  Eq.  (11): 

k  y  cos2  7]  cos  ( n  ,  z  )dA  ~  z  cos2 v  cos(n»y)dA] 


r  -  _ 

Si 


From  Fig.  C,  y  =  ±  -y-  cos  A  and  Z  =  0  ,  where  the  upper  sign  applies 
to  the  right  side  and  the  lower  sign  applies  to  the  left  side.  Since  y  is 
not  a  function  of  <£,  the  rolling- moment  coefficient  for  either  side  is 

ct  =  ±  -l-c;-a  (-cn) 


The  total  rolling- moment  coefficient  for  both  sides  is 

Ci  =  ~  2  7  A'  (CNl.ft  -  CNriBht) 


and  the  derivative  is 

r1  L  cos  A 


2  i 


(CNft.n 


“  CNo  ) 
Prigh  t  / 


For  either  side,  Eq.  (34)  gives 


-/  (  a  °°s‘  V  )  sin^A# 

v  <?/9  //3  =  0 

Integrating  gives 


CN0  = 


'P 

Thus, 


2  K  L  R 
3  S 


r 


A  sin  a  sin3  <£>  —  cos  a  cos  A  cos3  <b 


CNn  =  "CNn 

P  left  Pright 


and 


Cip  “ 


—  2  K  L  R  cos  A  sin  A 
3  Si 


£  sin  a  sin3  (f>  — 


cos  a  cos  A  cos 


*] 


(51) 


4>' 

tf>o  (or 


(52) 
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For  0  <  a  <  aQ 


C  t 


P 


S  i 


KL2R 


sin  a  cos  A  sin  A  sin  <f>' 


and  for  a0  <  a  <  ( tt  -  a0  ) 


Ci 


s  I 


(3  KL2R 


A  sin  A  J^sin  a  (sin3  cf> '  -  sin3  ^60) 

—  cos  a  cos  A  (cos3  tj>'  —  cos  0O)J 


For  the  hemicylinder  (<£'  =  7r  / 2 ) 


Ci 


Si 


/3  K  L2  R 


- —  cos  A  sin  A  sin  a  (  1  + 


a(: 


cos  a  cos 


•y/  1  —  s  in  2  A 


T) 


The  numerical  evaluation  of  Eq.  (55)  is  presented  in  Fig.  4d. 

2.1.4  Flat-Topped  Swept -Cylinder  Leading  Edge 


(53) 


(54) 


(55) 


The  geometry  of  the  flat- topped  swept- cylinder  leading  edge  is 
shown  in  Fig.  E. 


Fig.  E  Flat-Topped  Swept-Cylinder  Leading  Edge 


16 


AEDC-TDR-64-25 


The  directional  cosines,  pressure  coefficient,  and  elemental 
area  are  the  same  as  for  the  complete  cylindrical  leading  edge.  For 
0  <  a  <  77/2  ,  the  limits  of  integration  on  cf>  are  from  0  to  <f>' .  Then 
Eqs.  (34),  (39),  (45),  and  (52)  give 


Cn 


KLR 


(sin2  a  —  cos2  a  cos2  A )  (cos3  <f> '  —  1 )  ~  3  sin  a  (  cos  <j> '  —  1  ) 


+  2  sin  a  cos  a  cos  A  s 


in3  cf> '] 


(56) 


CA 


KLR 


cos  A 


(sin2  a  —  cos2  a  cos2  A)  sin  <j>'  +  3  cos2  a  cos2  A  sin  (f)' 


(57) 


CY 


/S  KLR 


*/3  KL2R 


—  2  sin  a  cos  a  cos  A  (cos  <f>'  —  1)J 
iin2  A  £  cos  a  cos  A  ( 3  sin  </> '  —  sin3  (/>')  —  sin  a  ( cos3  <f>'  —  1  )  J  (58) 


—  cos  A  sin  A  s in  a  sin3  '  —  cos  a  cos  A  (cos3  cf> 


)]  (59) 


For  the  hemicylinder  with  =  jt/2  , 


KLR 


(  •  2 
(  sin  a 


+  sin  a  cos  a  cos  A  + 


cos  a  cos 


(60) 


Ca 


KLR 


- —  cos  A 
3 


+  sin  a  cos  a  cos  A  +  cos  a  cos 


52  A 


(61) 


s 

KLR 


4 

3 


A  (< 


+  2 


cos  a  cos 


A) 


(62) 


r  ‘  s  i 

KL2  R  “ 


1 

3 


cos  A  sin  A  (sin  a  +  cos  a  cos  A) 


(63) 


The  pitching- moment  coefficient  and  yawing-moment  coefficient  de¬ 
rivative  are  given  by  Eqs.  (43)  and  (50),  respectively,  with  Cn  and 
Cy^g  determined  from  Eqs.  (56)  and  (58)  or  (60)  and  (62).  For  a  >  n/2  , 
the  equations  for  the  complete  leading  edge  also  apply  to  the  flat-topped 
leading  edge. 


The  characteristics  of  the  flat-topped  swept-cylinder  with  </>'  =  n/2 
are  given  in  Fig.  5. 
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2.1.5  Swept-Wedge  Wing 

A  planar  wing  having  sweepback  and  dihedral  is  analyzed  in  this 
section.  The  nomenclature  used  in  the  derivation  is  shown  in  Fig.  F. 


Fig.  F  Swept-Wedge  Wing 


The  sweepback  angle,  A ,  is  taken  as  positive  in  the  equations  for  the 
coefficients  of  the  entire  wing.  Where  the  sides  must  be  considered 
separately,  as  in  Eqs.  (65)  and  (66),  the  sign  convention  is  A  >  0  for 
the  right  side  and  A  <  0  for  the  left  side.  The  centerline  angle,  e ,  is 
always  positive  and  is  related  to  the  dihedral  angle,  r  ,  and  the  sweep- 
back  angle,  A ,  by 


€ 


(64) 


The  direction  cosines  of  the  inward  directed  unit  normal  vector  are 


cos  (n  ,  x ) 
cos  (n  ,  y  ) 
cos  ( n  ,  z ) 


—  sin  e _ 

yj  1  +  tan 2  A  sin2  f 

—  tan  A  sin  e 

\/  1  +  tan  2  A  sin  2  € 

_ ±  cos  € _ 

\/  1  4-  tan2  A  sin2  £ 


(65) 
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where  the  upper  sign  applies  to  the  upper  surface  and  the  lower  sign 
applies  to  the  lower  surface.  Then,  from  Eq.  (5) 

cos  B  sin  (e  ±  a)  +  sin  B  tan  A  sin  e  / 

cos  r]  -  - — — ■■■  -  . -  . - - -  '■ . : — - 

\J  1  +  tan2  A  sin2  e 

where  (e  +  a)  is  used  for  the  lower  surface  and  (<r  -  a)  is  used  for  the 
upper  surface.  At  /8  =  0  ,  the  upper  surface  becomes  shielded  from  the 
flow  at  a  =  e ,  and  the  lower  surface  is  shielded  at  a  =  n  -  e .  Equa¬ 
tion  (66)  is  valid  only  within  these  limits. 


Since  the  pressure  coefficient  is  constant  over  each  of  the  wedge 
surfaces  it  is  not  necessary  to  integrate  to  obtain  total  loads.  The 
force  and  moment  coefficients  are  most  easily  calculated  by  using  the 
projected  planform,  base,  and  side  areas.  The  total  planform  area, 

Ap  ,  is 

Ap  =  be  (67) 

The  base  area,  Ab  ,  for  either  half  of  wing  is 

Ab  =  bt  (68) 

The  side  area,  As  ,  for  either  half  of  wing  is 


In  the  following  derivations  it  is  assumed  that  Cp  =  0  on  the  base  of 
the  wing  and  that  the  wing  is  at  y3  =  0  .  Separate  equations  are  given  for 
the  lower  and  upper  halves,  and  the  total  wing  loads  may  be  obtained  by 
combining  the  two  halves.  No  graphical  results  are  presented  because 
of  the  simplicity  of  the  equations. 

2. 1.5.1  Normal-Force  Coefficient 


For  the  lower  half  of  the  wing  at  0  <  a  <  {  n  -  e)  the  pressure  coef¬ 
ficient  is 


C 


pL 


K  cos 


2 


7  = 


K  sin2  (f  +  a) 
1  +  tan2  A  sin2  e 


and  the  normal-force  coefficient  is 


(70) 


(71) 


at  (n  -  (?)  <  a  <  7T,  Cnl  =  o 

For  the  upper  half  of  the  wing  at  0  <  a  <  e ,  the  pressure  coefficient  is 


C 


Pu 


K  sin  2  (e  —  a) 

1  +  tan  2  A  sin  2  e 


(72) 
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and  the  normal-force  coefficient  is 

CNu  =  -  CPu  (73) 

At  €  <  a  <  7T,  CNu  =  0. 

2. 1.5.2  Axial-Force  Coefficient 


For  the  lower  half  of  the  wing  at  0  <  a  <  (n  -  e) 


CAl  -  CPL 

Ab 

S 

(74) 

For  the  upper  half  at  o  <  a  < 

€ 

CAn  =  Cpu 

Ab 

S 

(75) 

where  CPl  and  CPu  are  given  in  Eqs.  (70)  and  (72) 

2. 1.5.3  Pitching-Moment  Coefficient 


For  the  lower  half  of  the  wing  at  0  <  a  <  (n  -  e) 

C  mL  =  CnL  -fi  CAl  T  (~3  +  h) 
or  from  Eqs.  (71)  and  (74) 

c“l  -  [i--‘a-(-r+h)J  C6) 

For  the  upper  half  at  0  <  a  <  e 

c-u  -  -T2-  [t-  "  tan  f  (f  -  +  h)]  <77> 

If  the  leading  edge  is  tangent  to  the  wing  surface,  h  is  given  by 

h  =  R  cos  y  (^8) 

where  R  is  the  leading  edge  radius  and  y  is  given  by  Eq  (33b).  If  the 
leading  edge  is  not  tangent  to  the  wing,  Eq.  (78)  is  not  valid  and  h  must 
be  determined  from  the  leading  edge  geometry. 

2. 1.5.4  Side-Force  Coefficient  Derivative 

For  the  lower  half  of  the  wing  at  0  <  a  <  (n  ~  <r) 

Cyl  =  a  CPL  (79) 
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where,  by  definition 

ACpr  =  ( Cp  -  Cp  ) 

FL  \  ■‘left  Pright /] 

and  from  Eq.  (66) 


AC 


PL 


—  4  K  cos  /3  sin  (3  tan  A  sin  €  sin  (  €  +  a) 
1  +  tan  2  A  sin  2  € 


Then, 


and 


dP  )  P 

4KAS  tan  A  sin  f  sin  (e  +  a) 
S  1  +  tan 2  A  sin 2  e 

For  the  upper  half  at  0  <  a  <  e 

4K  As  tan  A  sin  e  sin  (f—a) 


(Cy/s)l  "  “ 


(CyP 


U 


1  +  tan  2  A  sin2 


(80) 


(81) 


(82) 


(83) 


2. 1.5. 5  Yawing-Moment  Coefficient  Derivative 

For  the  lower  half  of  the  wing  at0<a<(7r-e) 

(c»>l  -  -rhr  Cyl  +  -r(cw.  -  cAu(t)L]  (84) 

From  Eqs.  (74),  (79),  and  (80) 

(CAright  “  CAle£t)L  =  ~  CYl  “XT" 

and  (85) 

-fri1-  -dnr) 

For  the  upper  half  at  0  <  a  <  <r 

(S)n  =  foh  -TT  0  ~  ~t£t)  (86) 

2. 1.5. 6  Rolling-Moment  Coefficient  Derivative 

For  the  lower  half  of  the  wing  at  o  <  a  <  (n  -  e) 

(Ci>L  -  F[-CYL(F+h)+  -f- <87) 
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From  Eqs.  (71),  (79),  and  (80) 


and 


(%)L  “  (Cy/s)l  £  [  3  (  tan^  F  l)  h] 


For  the  upper  half  at  0  <  a  1  e 


(%)u  =  -(c^)D  i[t(i=Tr  -1)  -hJ 


(88) 


(89) 


2.2  BODY  COMPONENTS 

Body  components  typical  of  lifting  re-entry  configurations  are  the 
spherical  segment,  cone  frustum,  and  circular  cylinder.  The  loads 
on  these  components  are  computed  in  the  following  sections.  The 
method  of  combining  components  to  give  a  complete  configuration  is 
described  in  Section  2.  3. 

2.2.1  Spherical  Segment 


The  spherical  segment  is  a  basic  nose  for  bodies  of  revolution.  The 
nomenclature  used  in  the  derivation  is  shown  in  Fig.  G. 


Moment 


Reference  Point 


Fig.  G  Spherical  Segment 
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The  direction  cosines  of  the  inward  directed  unit  normal  vector  are 

cos  (n,  x)  =  —  cos  0  I 

cos  (n,  y)  =  —  sin  0  cos  cp  )  (90) 

cos  (n, z)  =  —  sin  0  sin  (f>  f 

Then,  from  Eq.  (5) 

cos  v  =  cos  /3  (cos  a  cos  0  +  sin  a  sin  0  sin  (f>)  +  sin  y8  sin  0  cos  <f>  (91) 

At  yS  =  0  the  surface  becomes  shielded  from  the  flow  along  a  curve  de¬ 
fined  by  4>  =  <f>o  ,  where 

<f>  =  —  sin  1  ( - - - (92) 

°  \  tan  a  tan  0  J 

This  equation  is  valid  only  for  0  >  {n/ 2  -  a)  ,  since  no  shielding  occurs 
for  6  <  {n/ 2  -  a)  .  The  elemental  surface  area  is 

dA  =  R2  sin  0  d  0  d  (fi  (93) 

It  is  assumed  that  Cp  =  0  on  the  base  at  all  angles  of  attack,  and  the 
coefficients  and  their  derivatives  are  evaluated  at  yS  =  0  . 


2. 2. 1.1  Normal-Force  Coefficient 


The  normal-force  coefficient  is  given  by 
2 

Cn  =  — ^ ff  cos2  7]  sin2  6  sin  if)  d  (f>  d  6  (94) 

s  A 

Because  of  the  limitations  on  the  shielding  equation,  the  evaluation  of 
Eq.  (94)  must  be  treated  as  three  separate  cases.  In  each  case,  the 
integration  is  taken  over  the  right  side  of  the  body  and  the  result  is 
multiplied  by  2. 


(I)  0  <  a  <  U/2  -  0b) 

2  0J,  77/2 

Cn  =  - — — . —  /  /  cos2  tj  sin2  0  sin  <f>  d  <f>  d  0 

S  0  —  77/ 2 

which  gives 

Cn  „  ^  2  =  ~ 7T  cos  a  s^n  a 

K  R  2 


(id 


( 77/2  -  0b )  -  a  -  n/Z 


Cn 


2  K  R2 

s 


77/2  —a  77/2 

/  /  cos  77  sin  0  sin  <f>  d  cf>  d  0 

_  0  —77/2 

0b  77/2 

+  f  f  cos  77  sin  0  sin  (f>  d  <f>  d  0 

77/2  —  a  (f>o  _ 


(95) 


(96) 


(97) 
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(III)  tt/2  <  a  <  (  tt/2  + 

2  @b  ^/2 

C]y  =  —  ^ —  /  f  cos2  7]  sin2  0  sin  <f>  d  ef)  d  0  (99) 

^  a  — n/2  <b 

T  O 

This  equation  reduces  to  Eq.  (98). 


To  make  the  spherical  segment  compatible  with  the  cone  frustum, 
let  #b  =  n/2  -  8.  Then  from  Eqs.  (96)  and  (98),  for  0  <  a  <  8, 

CN  IF  =  cos  a  sin  a  cos*  $  (100) 

and  for  8  <  a  <  (n  -  8) 


Equations  (100)  and  (101)  were  evaluated  for  a  =  0  to  180  deg  and 
8  =  0  to  70  deg,  and  the  results  are  presented  in  Figs.  6a  and  b. 


2.2. 1.2  Axial-Force  Coefficient 

The  axial-force  coefficient  is  given  by 
2 

CA  =  —  / /  cos2  77  cos  6  sin  6  d<f>  dO 

'  S  A 

Evaluating  Eq.  (102)  for  the  three  cases: 

(I).  0  <  a  <  (tt/2  -  6b) 

2  6 b  TT/2 

— R —  J  C  cos2  „  cos  Q  s jn  Q  deb  dd 
S  0  —tt/2 

tt  (  sin  a  sin4  6^  2  4  _  2 

~  l - - -  —  cos  a  cos  C7b  +  cos  a 


CA  = 


so 

Ca 


K  R 


(102) 


(103) 

(104) 
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(II)  ( 77/2  -  6h)  <  a  <  v/2 

CA  = 


77/2  —  a  77/2 

i  f  cos  tj  cos  8  sin  6  d0  d  6 

0  — 77/2 


8  b  77/2 

r  r  2 

+  J  J  cos  Tj  cos  0  sin  0  d  <^>  d  8 

n/i  —  a  <j)Q 


(105) 


then 


Ca 


KR 


cos  a  cos 


1  /  cos  8h  \  /  sii 

\  sin  a  )  \ 


a  sin4  8\, 


—  cos2  a  cos4  0k  +  cos2  a  )  |  —  +  sin  1  (  - - - 

/  L  2  \  tan  a  tan  (9b  /  J 

cos  a  cos  ^  0  2  ^  N  ./~7~2  ( 

+  - - -  v,  i  —  o  cos  ob )  V  sm  a  —  cos  8 b  r 


(106) 


(III)  77/2  <  a  <  (77/2  +  0b) 


The  equation  for  this  case  reduces  to  Eq.  (106).  Letting  (9b  =  77/2  -  8 
Eqs.  (104)  and  (106)  give  for  0  <  a  <  8. 


CA 


KR2 

and  for  8  <  a  <(n  -  8) 

q  ,  r 

CA  -  s  -  1 


77  /sin  a  cos  5  2  •4c*.  2  \ 

—  I - - -  —  cos  a  sm  o  +  cos  a  ) 


(107) 


KR2 


cos  a  cos 


' 1  /  sin  8  \ 

\  sin  a  / 


+  (-gial  g  cosll .  _  cos*  a  sih<  s  +  cosa 


°) 


+  sin 


_l  cos  a  sin  8  /i  o  *_2  c>\  ,/  2  •  2  ^ 

+  - -  i.  1  —  d  sm  6)  V  sm  a  —  sm  6 


in"*  /-iSS-S.)] 

\  tan  a  )  J 

(108) 


For  (77  —  8)  <  a  <  77,  Ca  =  0  . 

The  numerical  evaluation  of  Eqs.  (107)  and  (108)  is  presented  in  Fig.  6c. 


2.2. 1.3  Pitching-Moment  Coefficient 

The  resultant  force  acts  through  the  center  of  curvature,  and  the 
pitching  moment  about  the  reference  point  is  zero. 

2.2.1 .4  Side-F  orce  Coefficient  Derivative 

A  general  relation  for  the  side-force  coefficient  derivative  for  all 
bodies  of  revolution  can  be  obtained.  Consider  an  axisymmetric  body 
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which  is  pitched  through  an  angle  a'  relative  to  the  free -stream  veloc¬ 
ity  and  is  then  rolled  through  an  angle  0  about  the  X-axis  as  shown  in 
Fig.  H. 


Fig.  H  Force  Coefficients  on  Body  of  Revolution  at 
Combined  Angles  of  Attack  and  Sideslip 


The  force  acting  normal  to  the  X-axis  in  the  X,  -plane  is  defined,  in 
coefficient  form,  as  (Cn)«  =  0  and  remains  constant  as  the  body  is 
rolled.  Then 


Cy  =  —  (Cn)^  =  0  sin  0 


(109) 


By  resolving  the  velocity  vector  along  the  body  X-,  Y-,  and  Z-  axes,  it 
can  be  shown  that 


or 


tan  $  = 


sin  $  = 


tan  /3 


tan  /3 


A /  tan  2  /3  +  sin  2 


Then, 


Cy  =  -  (Cn) 


tan  /3 


—  0  yj  tan2  (3  +  sin2  a 

Differentiating  with  respect  to  /8  and  then  letting  /3  =  0  gives 

(Cy/»)/U.  = 


(110) 

(110a) 

(111) 


-c, 


(112) 
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Thus,  Cy^g  may  be  obtained  from  Eqs.  (100)  and  (101).  At  a  =  0  , 
is  determined  by  substituting  Eq.  (100)  into  Eq.  (112)  and  then  letting 
a  =  0  .  Thus, 

(C^L0  -  ~  gKRi2Sc-°-?--  <113> 

The  numerical  evaluation  of  Eqs.  (112)  and  (113)  is  presented  in 
Figs.  6d  and  e. 

2.2. 1.5  Yawing-  and  Rolling-Moment  Coefficient  Derivatives 

As  was  the  case  with  pitching  moment,  the  yawing-  and  rolling- 
moment  coefficients  and  their  derivatives  are  zero  at  the  moment 
reference  point. 

2.2.2  Flat-Topped  Spherical  Segment 

The  geometry  of  the  flat-topped  spherical  segment  is  shown  in 
Fig.  I. 


Fig.  I  Flat-Topped  Spherical  Segment 
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The  direction  cosines,  pressure  coefficient,  and  elemental  area 
are  the  same  as  for  the  complete  spherical  segment.  For  0  <  a  <  n/2 
the  limits  of  integration  on  are  from  0  to  n/2  and  the  limits  on  0  are 
from  0  to  0b  •  Then,  from  Eqs.  (94)  and  (102),  with  <9b  =  (77/2  -  8)  , 


Cn 


K  R 


+  sF  S  cos  8  (2  CQS *  S  _  ! 
2  \ 


+  —  cos  a  sm  a  cos 
2 


4  8 


.2  10  ■  2  2  o 

sm  a  —  -  sm  a  cos  o 

3 


)] 


(114) 


K  R 


j_  r  tt 
2 


^  cos  a  (1  —  sin  §)  +  ~  sin  a  cos 


4  8 


+  cos  a  sm  a 

where  8  is  in  radians. 


-  —  §  —  sin  8  cos  8  +  2  sin  8  cos  5^ 


(115) 


Equations  (112)  and  (113)  apply  only  to  complete  bodies  of  revolu¬ 
tion,  and  Cy^for  the  flat-topped  spherical  segment  must  be  obtained 
from  Eq.  (9)  which  gives 

Cy  =  -  —  / /  cos2  77  sin2  6  cos  9S  d  9S  d  <9  (116) 


For  the  right  side  of  the  segment 

Cl/3  = 
where 


v  t)  2  ®  b  n  /  *  / 

K  R  f  f  '  °  cos  7 


SJ  „  J 

0  0 


sin2  6  cos  </>  d  <f)  d  6 


(117) 


d  COS  Tj  \ 

dP  //3  =  o 


^  7/3  =  0 

2  cos  </>  sin  6  (cos  a  cos  6  4-  sin  a  sin  6  sin  <f>)  (118) 


Substituting  Eq.  (118)  in  Eq.  (117)  shows  that  Cy^g  is  the  same  for  both 
sides  of  the  segment.  Then,  multiplying  Eq.  (117)  by  2  and  performing 
the  integrations,  with  <9b  =  {n/2  -  8), 

-y  cos  a  cos4  8  +  sin  a(^  —  8  —  sin  8  cos  <5  —  ~  sin  8  cos3  «)]  (H9) 

The  moment  coefficients  and  their  derivatives  are  zero  for  the  indicated 
reference  point.  For  a  >  n/2  the  equations  for  the  complete  spherical 
segment  also  apply  to  the  flat-topped  spherical  segment.  The  charac¬ 
teristics  of  the  flat-topped  spherical  segment  are  presented  in  Fig.  7. 


2.2.3  Hemisphere 


Although  the  hemisphere  is  a  limiting  case  of  either  the  spherical 
wedge  (A  =  90  deg)  or  the  spherical  segment  {8  =  Odeg) ,  the  body  is  of 
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enough  interest  to  warrant  a  presentation  of  the  equations.  Letting 
8  =  0  in  Eqs.  (101),  (108),  and  (112)  gives,  for  0  <  a  <  n  , 


Cn 

s 

KR2 

II 

sin 

a 

(1  +  cos  a) 

(120) 

CA 

S 

KR2 

_  tt 

8 

(1 

+ 

2 

cos  a) 

(121) 

Y/3 

S 

KR2 

=  -  ■ 

77 

4 

(1 

+  cos  a) 

(122) 

The  moment  coefficients  and  their  derivatives  are  zero  for  a  moment 
reference  point  at  the  center  of  curvature.  The  characteristics  of  the 
hemisphere  are  presented  in  Figs.  2  and  6. 

2.2.4  Flat-Topped  Hemisphere 


For  0  <  a  <  tt/2  and  8  =  0,  Eqs.  (114),  (115),  and  (119)  give 


Cn 

s 

_  n 

(1 

,  9 

•  2 

cos  a  sin  a  +  sin  a) 

(123) 

KR2 

8 

Ca 

S 

_  77 

(1 

4-  9 

cos  a  sin  a  +  cos2  a) 

(124) 

KR2 

8 

Zi 

CY(3 

S 

77 

(  cos 

a  +  sin  a) 

(125) 

KR2 

4 

The  moment  coefficients  and  their  derivatives  are  zero  for  a  moment 
reference  point  at  the  center  of  curvature.  For  a  >  n/2  the  equations 
for  the  full  hemisphere  apply.  The  characteristics  of  the  flat-topped 
hemisphere  are  presented  in  Figs.  3  and  7. 

2.2.5  Cone  Frustum 


The  cone  frustum  is  frequently  used,  as  a  nose  or  flare  section  of 
lifting  bodies.  The  nomenclature  used  in  the  derivation  is  shown  in 
Fig-  -T 
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The  direction  cosines  of  the  inward  directed  unit  normal  vector  are 

cos  (  n  ,  x  )  =  —  sin  8 

cos  (n,  y)  =  -  cos  S  cos  ^  (  (126) 

cos  (  n  ,  z  )  =  —  cos  8  sin  cf> 

Then,  from  Eq.  (5) 

cos  7]  =  cos  /3  (cos  a  sin  8  +  sin  a  cos  8  sin  <$>)  +  sin  cos  8  cos  <f>  (127) 

At  )8  =  o  the  surface  becomes  shielded  from  the  flow  along  a  line  de¬ 
fined  by  <f>  =  </>o  ,  where 

"  -si"“  (-SH-)  (128) 


This  equation  is  valid  only  for  a  >  8  ,  since  there  is  no  shielding  of  the 
surface  for  a  <  8  .  The  elemental  surface  area  is 

dA  =  r  dVr  U29) 

sm  o 

It  is  assumed  that  Cp  =  0  on  the  flat  surfaces,  and  the  coefficients  and 
their  derivatives  are  evaluated  at  /3  =  0  . 


2.2.5. 1  Norma  1-Force  Coefficient 


The  normal-force  coefficient  is  given  by 

Rb 

Cn  =  — — — — 5—  f  f  r  cos2  77  sin  <b  dr  d  cb  (130) 

S  tan  8  0  Rn 

Since  <f>0  is  not  a  function  of  r  ,  the  first  integral  may  be  evaluated  to 
give 

p  KLRb  (  1  +  £)  f  2  ‘jljjl  /ioi\ 

b-N  —  - Ys -  cos  T]  sin  0  dcp  (131) 

where  ^  =  Rn/Rb . 

Because  of  the  limitations  on  the  shielding  equation,  the  integration  of 
Eq.  (131)  must  be  treated  as  two  separate  cases.  In  each  case,  the 
integration  is  taken  over  the  right  side  of  the  body  and  the  result  is  multi¬ 
plied  by  2. 

(I)  0  <  a  <  8 


Cn 


KLRb  (1  +  Q 
S 


n/2 

f  cos  rj  sin  <j>  d <f> 
~n/2 


which  gives 


Cn 


s _ 

KLRb  (1  +  O 


—  tt  cos  a  sin  a  sin  8  cos  8 


(132) 


(133) 
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(II)  8  <  a  <  (n  -  8) 


77/2 


which  gives 
Cn 


KLRb  (1  +0 


„  KLRb  (l  +  {)  r  2  .  ,  A  , 

CN  -  - ^ -  J,  cos  r}  sin  cp 

b  9o 


=  cos  a  sin  a  sin  8  cos  8  —  +  sin 


(134) 


\-Z-  +  sin"‘  (-Sani_)] 

|_  2  \  tan  a  J  J 

/-2_gjn!a,cos2g  +  si^gcos^N  V  sin"  a  -  sin2  8 
\  3  sin  a  cos  0  / 


(135) 


At  -  <5)^  a  <  7T ,  Cn  =  0.  Equations  (133)  and  (135)  were  evaluated  for 
a  =  0  to  180  deg  and  8  =  0  to  40  deg,  and  the  results  are  presented  in 
Fig.  8a. 

2.2. 5.2  Axial-Force  Coefficient 


The  axial-force  coefficient  is  given  by 

Rb 

Ca  =  -jr~  f  f  r  cos2  77  dr  d (j> 
s  <f>  Rn 

Integrating  over  r  , 

Ca  =  - -  tan  8  J  cos  77  d<p 

Z  o  (p 

Evaluating  Eq.  (137)  for  the  two  cases: 

(I)  0  <  a  <  8 

n  KLRb  (  1  +  ^)  , 


n/  2 

r 

—  77/2 


which  gives 

CA  - $-7 - ^  (2 

A  KLRb  (1  +  <f)  2 

(II)  8  <  a  <  (n  ~  8) 


cos  a  sin 


CA  -  KLRV(1t-€t  tan  8  f 


77/2 


<t>o 


which  gives 

CA 


(136) 


(137) 


tan  8 


KLRb  (1  +7) 


COS2  T)  d<f) 

l 

(138) 

f  sin2  a  cos  8 ) 

(139) 

cos2  17  d  (ft 

(140) 

h  sin2  a  cos2  8)  \  +  sin 

/  tan  8  \  ~1 

L  2 

V  tan  a  J J 

sin2  a  —  sin2  8  J 

(141) 
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At  (it  -  8)  <  a  <  n  ,  Ca  =  0  .  The  numerical  evaluation  of  Eqs.  (139)  and 
(141)  are  presented  in  Figs.  8b  and  c. 

2.2.5. 3  Pitch ing-Moment  Coefficient 


The  pitching- moment  coefficient  is  given  by 
Rb 


Cm  = 


K 


S  i 


,  ,  r  (Rb  -  r)  cos  T)  sin  <b  ,  ,  , 

j  j  -  __  dr  d cp 

cf>  Rn  tan  b 


Rb 


-  /,  /  r2  cos2  tj  sin  <£>  dr  d<f> 
^  Rn 


(142) 


Integrating  over  r, 
K 


Cm  = 


Si 


— 

/  2  _ 

3  \ 

3 

1 

/  r  Rb 

-  —) 

r 

tan2.  8 

\  2 

3  / 

~T” 

Substituting  Eq.  (131)  into  Eq.  (143)  gives 


Cm  —  C  ]\j 


Rb 


i  tan  8 


2  (1-D 

3  cos  2  8  (l-<f2) 


f  cos2  T]  sin  (f>  d  (f>  (143) 

*.  * 


(144) 


Since  the  limits  of  integration  on  <f>  do  not  enter  into  this  derivation, 
Eq.  (144)  is  valid  for  0  <  a  <  n  . 


2. 2. 5.4  Side-Force  Coefficient  Derivative 

The  side-force  coefficient  derivative  is  given  by  Eq.  (112), 

Cyg  =  -  (145) 

H  sm  a 

and  Cy^g  may  be  obtained  from  Eqs.  (133)  and  (135).  For  a  =  0  , 
substituting  Eq.  (133)  into  Eq.  (145)  gives 

(cY/s)a=0  t l rb  nrfT  =  -  17  sin  8  cos  s  (i46) 

The  numerical  evaluation  of  Eqs.  (145)  and  (146)  is  presented  in  Fig.  8d. 

2. 2. 5. 5  Yawing-Moment  Coefficient  Derivative 

It  is  obvious  from  Eq.  (145)  that  a  general  relation  for  all  bodies  of 
revolution  is 

^  cm 

^n{3  -  “  sin  a  (147) 

Substituting  Eqs.  (144)  and  (145)  in  Eq.  (147)  gives 
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Cn 


P 


=  CY/3 


Rb 


i  tan  8 


1  - 


3  cos 


( i  -  e ) 
( i  -  i2 ) 


(148) 


Like  Eq.  (144),  this  equation  is  valid  for  0  <  a 


2. 2. 5. 6  Rolling-Moment  Coefficient  Derivative 

The  resultant  force  acts  through  the  center  of  the  cone,  and  there 
is  no  rolling  moment  about  the  indicated  reference  point.  There¬ 
fore,  Cjt  =  0  . 

2.2.6  Flat-Topped  Cone  Frustum 


The  geometry  of  the  flat-topped  cone  frustum  is  shown  in  Fig.  K. 


i  =  Rn/Rb  t  Z 

Fig.  K  Flat-Topped  Cone  Frustum 


Y 


The  direction  cosines,  pressure  coefficient,  and  elemental  area 
are  the  same  as  for  the  complete  cone  frustum.  For  o  <  a  <  n/2,  the 
limits  of  integration  on  0  are  from  0  to  n/2  .  Then,  from  Eq.  (131)  and 
(137) 


Cn  - =  —  cos  a  sin  a  sin  8  cos  8  +  cos2  a  sin2  8  +  sin2  a  cos2  8  (149) 

1N  K  L  Rb  (  1  +  <f)  2  3 


and 


Ca  -  — . .  ^7 - =  tan  8  2  cos  a  sin  a  sin  8  cos  8 

KLRb  (1+  L 

+  ^cos2  a  sin2  8  +- 


sin  a  cos 


)] 


(150) 


Equations  (145)  and  (146)  apply  only  to  the  complete  cone  frustum,  and 
for  the  flat-topped  cone  frustum  must  be  obtained  from  Eq.  (9). 


33 


AEDC-TDR-64-25 


For  the  right  side  of  the  cone  frustum. 


and 


where 


Cy  =  - 

Cyp  =  - 


n/2  Rb 
K  r  r 

S  tan  8  o  Rn 

KLRb  (1  +  <f)  77/2 

2  S  fj 


cos  5  cos  cfj  (cos  a  sin 


2 

COS  T] 


cos  dr  d  <j> 


(151) 


cos  (f)  d  (f>  (152) 


8  +  sin  a  cos  8  sin  )  (153) 


Since  Cy^  is  the  same  for  both  sides  of  the  cone  frustum,  Eq.  (152)  is 
multiplied  by  2  and  integrated  to  give 


Cy/s 


_ S _ 

K  L  Rb  (  1  +  O 


—  cos  a  sin  5  cos  5 
2 


_2 

3 


sin  a 


(154) 


The  pitching- moment  coefficient  and  yawing-moment  coefficient 
derivative  are  given  by  Eqs.  (144)  and  (148),  respectively,  with  Cn 
and  Cy^g  determined  from  Eqs.  (149)  and  (154).  The  rolling- moment 
coefficient  derivative  is  zero.  For  a  >  n / 2  the  equations  for  the  com¬ 
plete  cone  frustum  apply  to  the  flat-topped  cone  frustum.  The  charac¬ 
teristics  of  the  flat-topped  cone  frustum  are  presented  in  Fig.  9. 

2.2.7  Circular  Cylinder 


The  circular  cylinder  is  a  special  case  of  the  cylindrical  leading 
edge  (A  =  90  deg)  and  the  cone  frustum  ( 8  =  0  and  £  =  1).  Letting 
A  =  90  deg  in  Eqs.  (38),  (42),  (43),  (49),  (50),  and  (55)  gives 


Cn 

s 

4 

.  2 

sin  a 

(155) 

KLR 

3 

CA 

=  0 

(156) 

cm 

=  Cn 

L 

2  i 

(157) 

CY  p 

S 

_4 

(158) 

KLR 

3 

sin  a. 

C”l8 

=  Cy/s 

L 

2Jt 

(159) 

C*R 

=  0 

(160) 

where  the  moment  coefficients  are  referenced  to  the  base  of  the  cyl¬ 
inder.  The  characteristics  of  the  circular  cylinder  are  presented  in 
Figs.  4  and  8.  The  Eqs.  (155)  through  (160)  also  apply  to  the  flat- 
topped  circular  cylinder. 
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2.3  COMPOSITE  CONFIGURATIONS 

The  Newtonian  analysis  is  based  on  the  local  flow  deflection  angle 
and  assumes  that  the  only  interference  between  components  is  that  due 
to  shielding.  Therefore,  the  configuration  being  analyzed  can  be  broken 
into  independent  elements  corresponding  to  the  components  described  in 
the  previous  sections.  From  the  equations  and  charts,  the  aerodynamic 
characteristics  of  each  component  can  be  determined  based  on  the  refer¬ 
ence  area  and  length  of  the  complete  configuration.  The  contributions  of 
the  individual  components  are  added  together  to  give  the  total  coefficients 
of  the  configuration.  When  this  method  is  used  to  obtain  the  character¬ 
istics  of  a  composite  configuration,  one  component  may  be  in  a  position 
where  it  shields  another  component  from  the  flow,  and  the  effects  of  this 
shielding  must  be  considered. 


The  summation  of  moment  coefficients  and  their  derivatives  re¬ 
quires  a  moment  transfer  from  each  component  reference  point  to  a 
common  reference  point.  Let  Xt  and  Zt  be  the  moment  transfer 
distances,  positive  as  shown  in  Fig.  L. 


Reference  Point  of  Complete  Configuration 


Fig.  L  Moment  Transfer  Lengths 


Then, 


c 

mcon  f  ig  ura  tion 

—  C  —Pat  XT 

mcomponent 

r 

+  C  A  , 

(161) 

C no 

Pc  onfiguration 

^component  ^  P 

xT 

1 

(162) 

,  ■ 

/-’configuration 

=  C‘A»p..„,  +  Cy/3 

ZT 

i 

(163) 

The  characteristics  of  a  typical  composite  configuration  are  evalu¬ 
ated  in  the  Appendix. 
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APPENDIX 

APPLICATION  OF  METHOD  TO  A  TYPICAL  DELTA  WING 


As  an  example  of  the  use  of  the  equations  and  charts  given  in  this 
report,  the  aerodynamic  characteristics  of  a  typical  delta  wing  were 
computed.  The  configuration  which  was  analyzed  is  shown  in  Fig.  10, 
and  the  lengths,  angles,  and  areas  used  in  the  calculations  are  given 
below: 

(I)  Spherical- wedge  nose  component 

R  =  0.0212  Ld  A  =  75  deg 

XT  =  -0.518  Ld 
Zt  =  0 

(II)  Swept -cylinder  leading  edge  component 

R  -  0.0212  LD  A  =  75  deg 

L  =  0.950  Ld  cf>'  =  74.5  deg 

Xt  =  0.40  Ld 
Zt  =  o 

(III)  Swept -wedge  wing  component 

c  =  0.940  Ld  A  =  75  deg  Ap/S  =  0.888 

t  =  0.0675  Ld  F  =  15  deg  Ab/S  =  0.064 

h  =  0.0205  LD  £  =  4.11  deg  As/S  =  0.119 

Xt  =  0.40  Ld 
Zt  =  o 

The  coefficients  were  based  on  the  planform  area  (S  =  0.267  Ld2)  ,  with 
the  mean  aerodynamic  chord  as  reference  length  ( i  =  0.667  Ld)  for  the 
pitching -moment  coefficient  and  the  span  as  reference  length  ( l  =  0.536  Lj) 
for  the  yawing-  and  rolling -moment  coefficients.  For  all  angles  of  attack 
it  was  assumed  that  the  base  of  the  delta  wing  contributed  no  axial  force, 
i.  e.,  Cpb  =  0  .  The  modified  form  of  the  pressure  coefficient  given  by 
Eq.  (2)  was  used  with  K  =  2  . 
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The  characteristics  of  the  nose  component  were  determined  from 
Fig.  2  and  the  moment  transfer  equations,  Eqs.  (161),  (162),  and  (163). 
Since  <f>'  was  not  equal  to  90  deg,  it  was  necessary  to  use  the  equations 
of  Section  2.  1.  3  and  the  moment -transfer  equations  to  evaluate  the 
leading- edge  component.  The  wing  component  was  evaluated  with  the 
equations  of  Section  2.  1.  5  and  the  moment-transfer  equations. 

The  aerodynamic  coefficients  of  each  component  and  of  the  complete 
delta  wing  are  presented  in  Fig.  11.  To  provide  a  comparison  of  the 
theory  with  hypersonic  experimental  results  for  this  delta  wing,  data 
obtained  in  the  50-in.  Mach  8  Gas  Dynamic  Wind  Tunnel,  Hypersonic  (B) 
of  the  von  Karman  Gas  Dynamics  Facility,  Arnold  Engineering  Develop¬ 
ment  Center, are  also  given.  The  data  were  obtained  at  a  Mach  number 
of  8.  1  and  free-stream  Reynolds  numbers  of  1.  3  to  5.  2  x  10®  based  on 
model  length.  The  theory  predicts  the  experimental  results  with  good 
accuracy  at  angles  of  attack  up  to  about  5  7  deg,  where  the  shock  wave 
becomes  detached.  Above  this  angle  of  attack,  the  theory  will  give 
better  agreement  with  the  experimental  data  if  K  =  Cp  =  1.83  is  used. 

m  y 
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Fig.  2  Aerodynamic  Characteristics  of  Spherical-Wedge  Noses 
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Fig.  3  Aerodynamic  Characteristics  of  Flat-Topped  Spherical-Wedge  Noses 
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Fig.  4  Aerodynamic  Characteristics  of  Swept-Cylinder  Leading  Edges  (<f>  =  jt/2) 
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Fig.  6  Aerodynamic  Characteristics  of  Spherical  Segments 
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c.  Axial  Force 
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d.  Side-Force  Derivative  (j 8  =  0),  5  =  0  to  45  deg 
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e.  Side-Force  Derivative  ( /S  =  0),  8  -  45  to  70  deg 
Fig.  7  Concluded 
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i.  Rolling-Moment  Derivative,  (3=0 
Fig.  11  Concluded 
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